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Abstract Some exact solutions and multi-mode invariant submanifolds were found for the 
Fermi-Pasta-Ulam (FPU) /3-model by Poggi and Ruffo in Phys. D 103 (1997) 251. In the 
present paper we demonstrate how results of such a type can be obtained for an arbitrary N- 
particle chain with periodic boundary conditions with the aid of our group-theoretical approach 
[Phys. D 117 (1998) 43] based on the concept of bushes of normal modes in mechanical systems 
with discrete symmetry. The integro-differential equation describing the FPU-a dynamics in 
the modal space is derived. The loss of stability of the bushes of modes for the FPU-a model, 
in particular, for the limiting case N ^ oo for the dynamical regime with displacement pattern 
having period twice the lattice spacing (7r-mode) is studied. Our results for the FPU-a chain 
are compared with those by Poggi and Ruffo for the FPU-/? chain. 
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1 Introduction 

Fermi, Pasta and Ulam (FPU) introduced in [|l| a simple nonlinear model for studying 
the problem of equipartition of energy among different degrees of freedom (linear normal 
modes) of an A^-particle mechanical system for sufficiently large A^. This model represents 
a monatomic chain with interactions only between neighboring atoms. It may be also 
considered as a chain of identical masses connected with identical nonlinear springs. The 
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force F{Ax) produced by each spring can be expressed as a power series in its deformation 
Ax: 

F{Ax) = -k-Ax + a- {Axf + (3 ■ {Axf + ... (1) 

The FPU-a chain corresponds to the case a 7^ 0, /5 = 0, while the FPU-/? chain corre- 
sponds to the case a = 0, /5 7^ 0, and in the both models all terms of higher order than 
those written in (|I]) are neglected. 

The FPU model played an important role in the development of some new concepts in 
the nonlinear dynamics of classical systems and helped to reveal a number of new nonlinear 
phenomena (see, for example, 0, |^, ^] and references cited therein). Let us specifically 



refer to the approximately recurrent behavior in time of the energy distribution among 
several of the first modes fl]], to the introduction of the concept of solitons @], to some 
important results concerning deterministic chaos 0, to breathing self- localized solitons 
i, etc. 

An interesting paper devoted to the FPU-/5 chain dynamics was published recently 
by Poggi and Ruffo 0. These authors revealed some exact solutions for the considered 
mechanical system and demonstrated "the existence of subsets of normal modes where 
energy remains trapped for suitable initial conditions" , which they called "subsets of I- 
type"[]. They also discussed the problem of stability of such solutions and appropriate 
multi-mode invariant submanifolds (with particular attention to the case of the thermo- 
dynamical limit 00). 

In connection with the above mentioned paper we want to emphasize the following 
essential points: 

• The analysis in |^ is based on the specific character of interaction between particles 
of the FPU-/? chain and, therefore, such an analysis must be done once again for every 
other type of monatomic chain. 

• Only one- and two-dimensional subsets of normal modes of "type I" were found for 
FPU-/5 in 0. 

In the present paper we would like to demonstrate a simple group-theoretical method, 
based on our previous papers [7-10], for finding all subsets of modes of "type I" for 
^We use another term for such objects - bushes of normal modes [7 - 10]. 
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monatomic chains with periodical boundary conditions. This method, in its general form, 
can be used for any mechanical system with any space symmetry group [0, |n|.0 It is 
based on the concept of "bushes of normal modes" in nonlinear physical systems with 
discrete symmetry, which we consider to be as fundamental as the concepts of solitons, 
dissipative structures, etc. in modern nonlinear science. Apparently, bushes of modes 
play an important role in many physical phenomena of current interest . For example, 

• pecuharities of dielectric and optical spectra brought about by the interactions be- 
tween different modes in crystals; 

• critical dynamics, behavior of structural and thermodynamical parameters near 
phase transition points; 

• phenomena induced by interaction between soft and hard modes, in particular, pecu- 
liarities of the temperature dependence of the Debye -Waller factor near phase transitions, 
which plays an important role in the interpretation of x-ray and neutron structural ex- 
periments, etc. 

The present paper is devoted to the symmetry- determined hushes of vibrational modes 
in the monatomic chains. Let us explain these basic notions on the qualitative level. The 
more exact definitions we give in the next section. 

Modes. We consider the symmetry group Go of a given mechanical system in equi- 
librium ("parent group") and its irreducible representations (irreps). The basis vectors 
of these irreps, which describe the atomic displacements patterns, are called "symmetry- 
adapted coordinates". In a particular case (the monatomic chain is precisely this case!), 
symmetry-adapted coordinates can be identical with normal coordinates. We mean by the 
term "mode" an arbitrary superposition (with time-dependent coefficients) of the above 
basis vectors of a given multidimensional irrep. All that was just said allows us to use for 
the monatomic chain the term "modes" (or "vibrational modes" ) as a synonym of both 
the symmetry-adapted and normal mode^. 

In the most part of the paper, we treat Gq as a cyclic group of pure translations 

^ In particular, it can be used for multiatomic chains. 

^ In more complicated cases when a given irrep enters several times into the full vibrational representa- 
tion of the considered mechanical system, we must differentiate between the terms "symmetry- adapted" 
and "normal" coordinates (modes). 
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(only in Sec. 2.6, we consider Gq as the dihedral group). All irreps of this group are 
one-dimensional, and each vibrational mode is simply a product of some time-dependent 
coefficient with the appropriate normal coordinate which can be introduced in the har- 
monic approximation by the conventional method. 

Bushes of modes. Normal modes are independent of each other in the harmonic ap- 
proximation, but the excitation from the initially excited mode (we call it "root" mode) 
can spread to a number of modes with zero amplitudes at the initial instant (we call them 
"secondary" modes), if anharmonic terms in the Hamiltonian of the considered mechanical 
system are taken into account. 

It is essential that, in general, not all modes turn out to be excited as a result of 
the excitation of the given root mode, but only a very certain their collection. This phe- 
nomenon is a consequence of the specific selection rules for the excitation transfer between 
normal modes of different symmetry[7]. The bush of normal modes is a superposition of 
all modes, associated with different irreps of the parent group Go, which are involved in 
the vibrational process as a result of exciting a given root mode. 

Every bush describes an exact dynamical regime whose number of degrees of freedom 
can be essentially less than that of the considered mechanical system.^ Note that one- 
dimensional bushes represent the similar nonlinear normal modes introduced by Rosenberg 
[15]. 

As was already stated, excitation of a primary ( "root" ) mode leads to the excitation 
of the bush as a single dynamical object. Amplitudes of modes belonging to the bush 
change in time and, as a rule, their evolution is not trivial. As a consequence, we can 
speak about a number of new types of excitations in systems with discrete symmetry. 

The following propositions were justified in the previous papers [7-10]: 

1. A certain subgroup G of the symmetry group Go corresponds to a given bush, and 
this bush can be excited by imposing the appropriate initial conditions with the above 
symmetry group G C Gq. 

2. Each mode belonging to the bush possesses its own symmetry group which is greater 

than or equal to the group G of the whole bush. 

Only in the special case, a given bush can possess the trivial symmetry and represent the general 
dynamical regime whose dimension coincides with that of the considered mechanical system. 
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3. In spite of evolving mode amplitudes, the complete collection of modes in the given 
bush is preserved in time and, in this sense, the bush can be considered as a geometrical 
object. 

4. The energy of the initial excitation is trapped in the bush, i.e. it cannot spread to 
the modes which do not belong to the bush, because of the symmetry restrictions. 

5. As an indivisible nonlinear object, the bush exists because of force interactions 
between the modes contained in it. 

6. Taking into account the concrete type of interactions between particles of the 
considered mechanical system can only reduce the dimension of the given bush. 

7. The extension of the bush can be realized as a result of the loss of its stability 
which is accompanied by spontaneous breaking of the bush symmetry (dynamical analog 
of phase transition). 

In the present paper, we discuss symmetry-determinedhnshes of vibrational modes, i.e. 
bushes whose existence is brought about by the symmetry-related causes only. In other 
words, we do not take into account any information about the type of interactions between 
the particles of the mechanical system. Nevertheless, it should be noted that taking into 
account such information can lead to some additional selection rules for excitation transfer 
between different modes and, therefore, to the reduction of the dimension of the symmetry- 
determined bushes (we discuss this problem in Sec. 2. 5). Moreover, every additionally 
considered symmetry can result in a similar effect (see Sec. 2. 6.). Thus, one can be sure 
that the excitation of a given mode cannot spread out of the symmetry-determined bush 
having that mode as root, but there may exist additional causes actually restricting the 
excitation to a mode subset of this symmetry-determined bush. 

Let us stress that the symmetry-determined bushes which are found in the present 
paper are valid for any of monatomic chains and, in some a sense, they can be applied to 
multiatomic chains as well (we explain this point in Sec. 2. 6). 

The notion of bushes of normal modes was introduced in [0]. In the last paper and in 
[8, 9] we investigated the main properties of bushes of modes and developed the general 
group-theoretical method for their finding. The detailed discussion of our approach and 



some important theorems are presented in |10|. Bushes of small dimensionality for many 



structures with point and space groups of crystallographic symmetry are found, classified 



and investigated in the above papers. Bushes of vibrational modes for the fullerene Ceo 
("buckyball" structure) are discussed in 



Since apphcation of group-theoretical methods for finding bushes of modes for mechan- 
ical systems, described by space groups with multidimensional irreducible representations, 
is very difficult, we developed a set of appropriate computer programs for treating such 
problems (previously similar programs were used for studying structural phase transitions 
in crystals 



All symmetry-determined "resonance subspaces" (see [T^, |T4|) and similar nonlinear 



normal modes were found in for all mechanical systems with any of the 230 space 
groups. Setting up this problem and its solution became possible because of using the 
concept of "irreducible" bushes of modes and due to the employment of the remarkable 
computer program ISOTROPY, created by Stokes and Hatch. This program realizes a 
great number of group-theoretical methods used in the theory of crystals. Its modified 
version, including the possibility to treat bushes of modes, is now available on the Internet 



as free software 16 



The first part of the present paper is devoted to the problem of existence of bushes 
of normal modes in nonlinear chains, and the second part is devoted to the discussion of 
their stability. Considering the existence of a bush of modes in nonlinear periodical chains 
with arbitrary interaction between their atoms (not necessarily of nearest-neighbor type), 
we use the general group-theoretical method described in |T0|. Discussing this problem 
we keep in mind not only the goal of obtaining some new results, but also the purpose of 
the exposition of our approach for studying dynamics of nonlinear systems with discrete 
symmetry using one-dimensional chains as the simplest case of such symmetry. Finally, 
we obtain for the FPU-a model a new integro-differential equation in the modal space 
for the case of the continuum limit (A^ oo). This equation is valid for all wave lengths 
(not only for long wavelengths as in the KdV-equation obtained for the FPU-a model by 
Zabusky and Kruskal |^). 

In the second part of the paper, we discuss the problem of bush stability using the 
FPU-a chain as a concrete model (in contrast to the FPU-/3 chain considered in [0). 
Different channels of the loss of bush stability for finite are discussed. The appearance of 
"satellites" forbidden by symmetry are revealed for bushes of normal modes for sufficiently 

6 



large values of A^. The simplest bush 5 [2a], corresponding to the displacement pattern 
with periodicity twice the periodicity of the chain in the equilibrium stateQ is studied for 
the limit N oo. 



2 Bushes of normal modes for nonlinear chains 



We describe here the general method for finding bushes of vibrational modes ||T^ using 
the one-dimensional FPU model with periodical boundary conditions as an illustrative 
example. The starting point of our approach is the expression of the set of atomic dis- 
placements as a sum of contributions from different irreducible representations (irreps) of 
the symmetry group of the considered mechanical system in equilibrium [see (|^) below]. 
Therefore, we begin with the consideration of irreps and their basis vectors for the case 
of the monatomic chain. 



2.1 Vibrations of monatomic chain in terms of irreducible rep- 
resentations of its symmetry group 

Let us consider the many-particle Hamiltonian system moving near the single equilibrium 
state characterized by the parent group Go of point or space symmetry and introduce 
the A'^-dimensional vector X(t) determining values of all of its degrees of freedom at a 
given instant t : 

X(t) = {xi(t),a;2(t),...,xjv(t)}. (2) 

We can write the vector X(t) as the superposition of basis vectors (p^/^"' of all those irre- 
ducible representations Tj of the group Gq which enter into the mechanical (reducible) 
representation of the considered system: 

x(t) = E/^i''^W-v'? = EA,(t). (3) 

jk j 

Here the index j corresponds to the irrep Tj while the index k corresponds to its different 
basis vectors because, in the general case, Tj can be a multidimensional representation. 
The iV-dimensional vector Aj(t) is the contribution to X(t) from the irrep Tj. 

^ This bush consists of only one zone boundary mode or 7r-mode, in another terminology. 
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For an A^-particle nonlinear chain, X()f:) is the complete set of displacements Xi{t){i = 
1,2,..., A^) of each particle (mass point) from its equilibrium position]^. [In accordance 
with the periodical boundary condition, we assume X7v+i(i) = xi(t)]. 

The translational group T can be chosen as a simplest variant of the parent symmetry 
group Gq. Nevertheless, the monatomic chain is invariant under the action of inversion 
with respect to its center and, therefore, the full symmetry in this case is the dihedral group 
D. Certainly, this group is more relevant for the symmetry-based analysis of nonlinear 
vibrations in monatomic chains. On the other hand, the inversion can be absent for some 
multiatomic chains, and the full symmetry for such cases is represented by the group T. 
Since some important results can be obtained for any one-dimensional chain with the 
aid of the group T only, we start our consideration with the more general and simple 
case Go = T. The influence of the different choice of the parent group [Gq = T or 
Go = D) on the group-theoretical analysis of nonlinear vibrations of monatomic chains 
will be considered in Sec. 2.6. 

Thus, we suppose that the equilibrium state of our monatomic chain is described by 
the parent group Gq = T which, obviously, happens to be the cyclic group of order 
with generator a: 

Go = {E, a, a^ a^ . . . , a^"^}, = E. (4) 

{E is the identity element of the group Go). Here d is the operator translating all particles 
by the space period a of the considered chain in equilibrium. This operator generates the 
cyclic permutation of all particles and, therefore, the permutation of the corresponding 
displacements Xi{t): 

d{Xi{t),X2it), . . . XN-lit),XNit)} = {xNit),Xi{t),X2it), . . . XN-l{t)} . (5) 



As is well known [^, the cyclic group Go of order A^ has A^ irreps and they are 
all one- dimensional!]. Moreover, the 1x1 matrices of its generator a for different irreps 
Ffc {k = 1,2,..., A^) are simply the A^th-degree roots of 1. As an example, we give all matrix 
irreps of the cyclic group Go for A^ = 12 in Table l|. (Since all are one-dimensional 
irreps. Table 1 coincides with the table of their characters). 
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We consider longitudinal vibrations of the chain. 

^In contrast to many point and space groups which can possess multidimensional irreps. 

®Note, that in this table and in some points of the text we neglect the hat over symbol a for simplifi- 



Let us construct the basis vectors of the irreps of the group Go in (4) which describe 
the different modes of vibration of our nonhnear chain. The basis vectors of the irreps are 
usually obtained by the method of the projection operator |l7|, but in our case it is simpler 



to make use of the "direct" method based on the definition of the group representation 
(see [|lOl)- Indeed, according to this definition 

a$ = M(a)$, (6) 

where M{a) is the matrix corresponding in the irrep to an element d of the group Gq, 
and $ is the complete set of basis vectors c^^ of this representation which are transformed 
under the action of the operator d [in we use tilde as a symbol of matrix transposition] . 

The matrix Mk{a) of the generator d of the cyclic group Go for the irrep {k = 
0,1,2, . . . , N — 1) is equal to 

Mk{a)=-f\ (7) 

where 7 = e^^^^^ and, therefore, 7^ = 1. 

The desired basis vector cpf, of the irrep F^ can be written in the form 

^, = {x['\xi'\...,x^^^), (8) 

where x\''^ {i = 1, . . . , N) are unknown displacements of the particles. From (^) and (|^, 
we get the following equation which can be used for obtaining the basis vector (p;^ : 

a(^, = 7Vfc, A; = 0,l,...,iV-l. (9) 

Equation (|^) leads us to a chain of simple algebraic equations for the components of the 
vector (fij^ in (H). Solving these equations, we obtain 

ifij, = (x, 7^ ■ X, f ^ ■x,f''-x,..., ■ x), (10) 

where x is an arbitrary constant determining the displacement of the first particle of the 
chain. (Hereafter, we denote complex conjugation by a bar over the appropriate value, 
for example, 7 = e~^'^*/^). The presence of only one arbitrary constant in the general 
form of each basis vector means that every irrep F^ of the group Go enters one and 
cation of the notation. 
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only one time into the decomposition of the mechanical representation of the monatomic 
chain. 

The above conclusion applies to the monatomic chain, but it is not valid for arbi- 
trary physical systems: a given irrep may enter several times or not enter at all into the 
appropriate mechanical representation. The method of finding the basis vectors for mul- 
tidimensional irreps and complex mechanical structures turns out to be very complicated 
and it is useful to use in this case the computer programs [l^, [l^ . 



Normalizing the vector c^^^ from Eq. (|ToD, we obtain the final form of the basis vectors 
of the irreps Tk of the group Go = T: 

^k^<p (r.) = -^(1, f , f\ f\ • • • , 7^^-')') • (11) 



These basis vectors provide us so-called symmetry- adapted coordinates. According to the 
Wigner theorem, they are normal coordinates associated with F^, if only one copy of 
this irrep (multidimensional, in general case) is contained in the decomposition of the 
mechanical representation into its irreducible parts^. 

As was shown above, the two notions — symmetry-adapted coordinates and normal 
coordinates — are identical when we characterize the monatomic chain by the translation 
group (Go = T). We will refer to the modes determined by the above coordinates as 
normal modes or vibrational modesQ. 

Each basis vector ^p^^ determines a certain pattern of displacements of all atoms of 
the chain and, in some cases, it is more convenient to consider real modes 'j/'fc instead of 
the complex modes f^^. 

Note that two basis vectors, for which the sum of their indices is equal to A^, are 
complex conjugates of each other: 



^N-k = Vk^ A; = 0,1,2, 



A^ 
~2 



(12) 



® Let us remember that in the opposite case, according to the Wigner theorem, additional hnear 

transformation is required to diagonahze the force constant matrix. In our general approach we use 

symmetry-adapted coordinates since they can be obtained with the aid of the group-theoretical methods 

alone (without any information about concrete interactions in the mechanical system). 

Modes differ from generalized coordinates y^, by the appropriate time-dependent functions ^k{t) 

in front of the vectors </3j, as in Eq. (^^. In further consideration, we will sometimes call by the term 

"modes" not only /ife(i) • y^., but also the vectors tp^. and even the functions fikit)- 
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(Hereafter, y denotes the integer part of y.) Indeed, 

^ ^^ ^XN-k) \ 1 /I --k ^-2k \ ^k ^2k \ — - 



/27rk \ /2nk ^ 

sm —--n + cos — r— 



n=l,2,...,iV . (13) 



since 7^ = 1 and 7"^^ = 7^^. 

Using the appropriate hnear transformation in the space of these mutually conjugated 
vectors, we obtain the real normal modes in the form given by Poggi and Ruffo : 

Here n is the number of the particle, and k {k = 0,1,2, . . . , N — 1) is the number of the 
mode. We call the real vectors i/))^ and V'Ar-fc "conjugate modes" similar to the term for 
pairs of complex modes ipf. and 'PM-k- 

Returning to the basic Eq. (§), we can rewrite it in the following form, 

N-l N-1 

x(t) = EM^)v'fc= E^fcWfc- (14) 

A:=0 fc=0 

Thus, the pattern of atomic displacements X(t) for the chain in the vibrational regime 
at any instant t can be represented as a superposition of complex ((/J^r,) or real (i/'fc) basis 
vectors with the coefficients ^k{t) or Vk{t) which depend on the time t. 

Let us emphasize that i/^jt from Eq. ([T3|) are normal modes of the linear chain (they 
can be obtained in the harmonic approximation), but we use them in ([T^ ) as the basis of 
the configuration space for studying dynamics of the nonlinear chain. 

2.2 Bushes of normal modes 

We continue to adapt the general method [7-11] of studying nonlinear dynamics of me- 
chanical systems with arbitrary discrete symmetry for the case of nonlinear chains. Every 
basis vector (normal mode) from Eqs. (pi]) or (0) determines a specific set of atomic 



displacements, and this pattern is characterized by a certain symmetry group Gk which is 
a subgroup of the group Gq of the chain in equilibrium {G^ C Gq). Let us excite a dynam- 
ical regime of the considered system using the set of atomic displacements corresponding 
to a given root mode V'fc ^is initial conditions at the instant t = 0: 

Xi\t=o = Aipki; Xi\t=o = ^ {i = l,2,...,N). (15) 
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Here ipki are components of the vector -0;, = {ipki, 4'k2, • • • , i^kN) and A is a fixed constant. 

The symmetry group of the mechanical system at an arbitrary later moment t cannot 
be lower than the group at the initial moment determined by the vector Indeed, 
the symmetry of the system at some instant t > to is determined by the intersection of 
the symmetry of its Hamiltonian and the symmetry of the initial conditions. On the other 
hand, Go can be considered as a symmetry group of the Hamiltonian, as well as of the 
equilibrium state, if this state is nondegenerate]^. Thus, every symmetry element g E Gk 
does not change the dynamical state X(t) for t > to {a. more detailed consideration of this 
problem can be found in [7,8]). 

Now we can pose the question: what modes Vfc from Eq. (]T^) can contribute to the 



considered dynamical regime X(t) corresponding to the group C G for t > to? Ob- 
viously, because of nonlinear interactions between different modes, the root mode 
brings about the appearance of a number of "secondary" modes in the considered dynam- 
ical regime, and we can obtain all these modes by demanding that the symmetry X(t) is 
equal to Gk (at least, it must not be lower than G^!): 

G,X(t)=X(t). (16) 

This equation means that X(t) is an invariant vector with respect to all elements of the 

group Gfc [5 As a consequence, only those modes ipj contribute to X(t) whose symmetry 

group Gj is higher than or equal to Gk of the initial configuration of the mechanical system 

in hand. The complete set of these modes, i.e., the root mode and all secondary modes 

corresponding to it, forms the hush as the geometrical object. 

Note that each subgroup G' C Go must be tried in finding bushes, but in general, 

some subgroups do not generate any vibrational hushes (see examples in Sec. 2.6). 

Let us consider the symmetry group Gk of different modes ily^ in more detail. We can 

We state this hypothesis at the beginning of Sec. 2.1. For the degenerate equihbrium state, i.e. when 

there exist several equivalent minimums of the potential energy, the symmetry group of the Hamiltonian 

is higher than that of the each of the appropriate equilibrium states. 

Rigorously speaking, we must distinguish between the symmetry elements g C Gk, which act on 

three-dimensional vectors of the Euclidean space, and the operators g, corresponding to them, which act 

on iV-dimensional vectors X(t) (see [10]). The complete set of the operators g associated with all g E Gk 

forms the group Gk- 
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find the group Gk C Go of a given mode by selecting all elements g G Gq under which 
this mode is invariant: 

9'il^k = '^k- (17) 

But there exists a more convenient method. Indeed, X(t) was written in @ as a sum of 
contributions Aj from different irreps Tj of the dimension rij. It is easy to show that the 
following consequence of Eq. (^) takes place for each individual irrep r^: 

G,A, = A,. (18) 

According to the definition (^ of the representation of the group Go we can act on the 
set of basis vectors i/?^"''' by the matrix M{g), instead of acting on them by the operator 
g {g & Go). This is a trivial procedure for the cyclic group Go, since all of its irreps are 
one-dimensional and, therefore, these matrices are scalar values.[^ 

On the other hand, Eq. (0) can be true if and only if the 1x1 matrix M{g), cor- 
responding to the element g E Go, is equal to the unit matrix. Therefore, all elements 
g E Go whose matrices of the irrep Tk satisfy the condition 

Mk{g) = 1 (19) 

do not change the mode t/?;, and, obviously, the complete set of such elements forms the 
subgroup Gk of the group G^ C Go. 

Let us consider the case = 12 as an illustrative example. All irreps of the corre- 
sponding cyclic group Go are given in Table 1, as well as the symmetry groups Gk of 
modes (pj. associated with these irreps. As was explained above, only those elements 
g E Go belong to the group Gk whose matrices in the irrep Tk are equal to 1. In Table 1, 
we define the subgroups Gk C Go by putting their generators in square brackets. Here- 
after, we write the element as {p ■ a): such notation is often used for cyclic groups. 
Note, that according to the well-known Lagrange theorem, the order (M) of any subgroup 
is a divisor of the order (A^) of this group, and that every subgroup of the cyclic group is 
also a cyclic group which can be determined by the single generator. Irreps rjv_fc and Tk 

are mutual conjugate representations: T^-k = Tk- It is easy to notice that all subgroups 
^•^In the general case, the appropriate procedure can be reahzed usmg the notion of "invariant" vectors 
of the irreps of the group Gq. 
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Gk of the original cyclic group Go with the order = 12 are presented in Table 1, and 
that the subgroups corresponding to the conjugate irreps and T^_k are identical (more 
exactly, subgroups corresponding to the conjugate modes and (p^^k associated with 
these irreps). Below we list subgroups Gk of the group Go with = 12, listing in braces 
all their elements (see footnote 8 on page |^), the order (M) and symbols of irreps whose 
basis vectors (modes) possess such a symmetry (G^). 

1. GoH = {E, a, 2a, 3a, ... , 11a}, Mq = N = 12. 

Only one mode (p^ = -^(1, 1, 1, • • • , 1), corresponding to the irrep To, possesses this 
symmetry group. The mode cp^ plays a special role: it corresponds to the movement 
of the chain as a whole. In the next sections, this mode is not taken into account 
since we consider only vibrational modes. 

2. Gi [2a] = {E, 2a, 4a, 6a, 8a, 10a}, Mi = 6. 

It can be found from Table 1 that only two modes </?o and </?g satisfy such symmetry, 
because the unit matrix (1) corresponds to the generator [2a] only for the irreps To 
and Fg. We write this fact as follows 

Gi[2a] :ro,r6. 

The symmetry group of the basis vector cpg is equal to Gi[2a], while that of the 
basis vector is higher (it is equal to Go [a]). Therefore, cpg is the root mode, and 
(Pq is a secondary mode corresponding to it. 

Proceeding in a similar way, we obtain: 

3. G2[3a] = {E,3a,6a,9a}, Ms = 4 : 10,14,18. 

4. G3[4a] = {E,4a,8a}, M3 = 3 : 10,13,16,19. 

5. G46a] = {E,6a}, M4 = 2 : 10,12,14,16,18,110. 

6. G5[12a] = {E}, M5 = l: for j = 0, 1, 2, . . . , 11. 

Now we can write all possible bushes for the chain with = 12 particles using Eq. ([T^) 
and the above results. Each bush corresponds to the appropriate subgroup Gk of the 
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parent group Gq, and it embraces all modes with symmetry equal to (for root modes) or 
higher than (for secondary modes) G^- 



Gi[2a] : 


x«(t) 




G2[3a] : 


x(^)(t) 


= l^o{t)^o + + f^8{t)<p8, (20) 


GsM ■■ 


x(3)(t) 


= f^o{t)^0 + /^3 {t) ^3 + (t) "Ps + f^9{t)y^9, 


Gi[6a] : 


x(^)(t) 


= f^o{t)^0 + f^2{t)V2 + ^4(t)¥'4 + f^6{t)^6 + f^8{t)^8 + f^lo{t)^io, 


G^lUa = E] : 


x(5)(t) 


11 

= ^^fij{t)(pj, (the trivial bush) . 



j=0 



Thus there exist four nontrivial vibrational bushes (besides the trivial bush with 
6*5 [12a] = E) for the nonlinear chain with = 12 particles. The modes (basis vec- 
tors of appropriate irreps) cp^. determine the pattern of atomic displacements, and the 
time- dependent coefficients fikit) characterize the evolution of these displacements in 
time. 

It easy to find from ( pOD that every bush with symmetry group Gk contains the con- 
tributions from modes with this symmetry Gk and from all other modes whose symmetry 
is higher than Gk- This is the geometrical cause of the spreading of the initial excitation 
from the root mode to other (secondary) modes with Gk equal to or higher than that of 
the root mode. The dynamical reason for such a phenomenon will be considered below in 
Sec. 2.3. 

Returning to the general case, we obtain that the numbers and the types of possible 
bushes for any nonlinear chain with particles depend essentially on divisibility prop- 
erties of the integer A^. Indeed, the total number of bushes for an A^— particle chain is 
equal to the number of divisors of A^. The number of modes contained in the concrete 
bush corresponding to the divisor M (this divisor determines the order of the appropriate 
cyclic group Gk) is equal to N/M.f^ In particular, for prime number A^ there exists only 

one (trivial) bush with Gi = {E} which contains all modes ipj {j = 0, 1,2, A^ — 1).[3 

^''Note, that /ife(i) associated with different bushes are not related. Indeed, we may introduce the 

additional index i for these coefficients ii'j^\t), which corresponds to the number of the bush, but this 

index is omitted for simplification of the notation. 

-•^^ Moreover, it is easy to verify that the order M of the group Gj[ma] equals N /m and, therefore, the 

number of modes of the bush with such symmetry is simply equal to m. 

In this case of trivial symmetry, the bush describes the full dynamics of the considered mechanical 
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Therefore, the bush structure for the chain with = 12 differs in principle from that for 
N =lloT N = 13. 

In conclusion, let us note that all possible symmetry-determined bushes of normal 
modes can be found for arbitrary mechanical system with discrete symmetry by the regular 
group-theoretical methods described in detail in our previous papers (see |10| and [p]). 
We will not consider these methods in the present paper, but we want to extract their 
main points. 

A certain subgroup G of the parent group Gq corresponds to a given bush B[G]. This 
subgroup must be fixed for singling out the complete collection of the normal modes 
contained in fifC]-!^ Then the structure of the given bush, i.e. the list of all modes 
belonging to -BIG], can be found as a result of solving the following systems of linear 
algebraic equations for each of the irreps Tj of the group Gq: 



(r, i G) C, = C,. (21) 

Here Tj | G is the restriction of the irrep Tj of the parent group Go to the subgroup 
G of this group. In other words, we must solve the equations 



for all matrices Mj{gi) corresponding to the generators Qi of the subgroup G in the irrep Tj. 
The invariant vectors Cj from (^) determine the set of coefficients in the superposition 
of the basis vectors of the given irrep Tj which is compatible with the symmetry group 
G of the given bush i?[G]. In particular, the invariant vectors can turn out to be zero 
vectors for some irreps Tj, and this means that such irreps do not contribute to the given 
bush. 

For obtaining the bushes of vibrational modesQ, in addition to the above results, 

we must find the explicit form of the basis vectors constructed from the set of atomic 

system without any simplification. 

In particular, the subgroup G can be fixed by indication of the initially excited root mode, since the 

definite group is associated with each normal mode. 

There exist bushes of modes of very different physical nature. 



16 



displacements, for all irreps which contribute to the given bush. (This can be done for 
example, with the aid of the method of projection operators). 

As we already noted, the above procedures are simple for the monatomic chain, but 
special group-theoretical computer programs are needed for the general case (see references 
m and PI). 

If the complete list of modes of a given bush is known, the root mode can be frequently 
defined as that with the lowest symmetry of all other bush modes. Nevertheless, there 
exist cases where the root mode of the bush can be chosen in different ways and there 
exist cases where the symmetry group of the bush is determined by the intersection of 
the symmetry groups of the modes associated with several different irreps of the parent 
group. In the last case, we speak about a bush with several root modes (see and, 
especially, and |^). In general, the set of root modes must uniquely determine the 
symmetry group of the whole bush. The excitation of this set of root modes leads by 
necessity to the excitation of all other modes of the bush. 



2.3 Bushes of modes as dynamical objects 

Considering bushes of modes in the previous section, we made use of geometrical (group- 
theoretical) methods only. Now let us consider the concrete mechanical model — the FPU- 
a chain described by the Hamiltonian 

1 ^ 1 ^ a ^ 

H = 7:J2pI + 7;T. i^n+i - Xnf + - E ' x^f • (22) 

^ n=l ^ n=l n=l 

Here x„(t) is the displacement of nth particle from its equilibrium state at the instant t, 
and Pn{t) is the corresponding momentum. 

As in the original paper by Fermi-Pasta-Ulam [|T|, the nonlinearity of the chain is 
assumed to be weak. It depends on two different factors — on the coefficient a in (|22| ) 
and on the energy of the initial excitation of the system. We can remove a from ( ^2] ) by 
means of a scaling transformation of x„(t), i.e., the coefficient a can be made equal to 1. 
Nevertheless, sometimes it is convenient to write a {\a\ -C 1) in Eq. (^2]) explicitly to 
stress the weakness of the nonlinearity of the chain. 

As was already stated, the basis vectors cpj (0) of irreps of the group Go are (linear) 
normal modes of the FPU chain in the harmonic approximation (a = 0), and we can use 
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them in (|T^ as a basis for studying nonlinear dynamics of the system with a 7^ 0. 
We consider the equations of motion in the form, 

dU 



dXr 



(n = l,2, 



,iV) 



(23) 



where U{xi, . . . ,xn) is the potential energy of the FPU chain. This system of differential 
equations becomes linear in the harmonic approximation (a = 0), and reduces to the set 
of independent oscillators 



i,,+uj]fij = 0, (j=0,l,2,...,iV-l) 



with frequencies 



2sm(- 



(24) 



(25) 



as a result of the transition from the initial variables x„(t) to the normal modes in the 
complex (|ll|) or real ([T3|) form. From this point, we call the time-dependent coefficients 
fij{t) [or z/j(t)] from Eq. (|I4D by the term, "normal modes", even though this term cor- 
responds really to fij{t)(pj. This is convenient because (pj are constant vectors and fij(t) 
are our new dynamical variables. Substituting ( [14|) into equations (|2^) and solving them 
with respect to variables /ij(t), we obtain the following nonlinear differential equations 
(see details of such a transition in the general case in the Appendix of Ref. ||10||): 



N-l 
k=0 



27r \ /27r 

2 sm — (1 — k) ] — sm — ? 



(26) 



7 



2ai 



N 



j = 0,l,2,...,N-l. 



Hereafter, we use the following convention about indices of normal modes: these indices 
must belong to the interval [0, — 1] and, therefore, the value of the mode index of 
the type ±j ± k {j,k = 0, 1, . . . , A^ — 1) must be reduced to this interval by adding or 
subtracting mN, where m is an integer number. 

The transition from initial variables x„(t) to normal modes i^jit) in the real form (|T3|) 
leads to the following dynamical equations for the FPU-a model: 



a 



N-l 



N ^0 



2n. 



2n 



[vj+k + v-j-k) ( 2 sin — + sin —j 



(27) 



27r, 2tx 
+ [^j-k - ^-i+k) ( 2 sin —A; - sin —j 



N 



j=0,l,2,...,Ar-l. 
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Now we can consider the dynamical cause of the existence of bushes of normal modesQ. 
First of all, let us remember that the mode /io corresponds to the movement of the FPU 
chain as a whole object and, therefore, to the movement of its center of mass. Vibrations 
associated with all other modes (j = 1,2,...,A^ — 1) cannot move the center of 

mass because of the specific type of displacement patterns described by basis vectors ipj. 
Therefore, the excitation of any mode yUj(t) with j ^ cannot lead to the excitation of 
the mode /zq, as a consequence of the conservation of momentum. This general conclusion 
can be easy verified using the dynamical equations (^61) . Indeed, we obtain the following 
equation of motion for the mode fiQ-. 

/io + 0/io = . 

It is also easy to verify that mode fio{t) cannot excite any other mode, and, therefore, we 
can forget about fio{t) when discussing dynamical equations. 

On the other hand, this mode enters into each bush of the nonlinear chain (see, for 
example, Eqs. (^)). What does it mean? We obtain bushes, as geometrical objects, from 
the idea of invariance of vibration patterns with respect to the subgroups of the cyclic 
group Go only. But there can exist further causes for reducing of the number of modes in 
a given bush because of some additional symmetry properties of the concrete dynamical 
system and because of its specific character (we discuss this problem in Sees. 2.5 and 2.6). 
In the above case of the mode fio{t), such cause is provided by the law of conservation of 
momentum which is the consequence of homogeneity of the space. 

Let us assume that only one mode fJ^Q{t) is excited in the FPU-a chain with = 12 
particles. In other words, we suppose that all modes except mode f^eit), are equal to zero 
identically (we call them "sleeping modes"): 

fij{t)=0, forjVe. (28) 

Can such a dynamical regime exist in our nonlinear chain? Using the above assumption, 
it is easy to obtain from Eqs. (^61) , 

jle + 4/i6 = 0, (29) 

llj + ujjfij = 0, for all j 7^ 6. 
^^The following analysis will be based on studying equations (|2^ ) for complex modes l^j{t), but the 
same results can be obtained using equations ( p7| ) for real modes Vj (t) as well. 
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Therefore, the hypothesis that ^eit) ^ while fij{t) = for j 7^ 6 is self-consistent, and 
we have a pure harmonic regime HQ{t) = Acos(-\/2 t + S) for the FPU-a chain {A,S are 
arbitrary constants). It corresponds to the one-dimensional bush {^Pq} with symmetry 
group Gi[2a]g [see Eqs. (|0D]. 

Let us consider the bush with symmetry [3a] corresponding to the second equation 
of (pop. We must now assume that fi^^ 0, /is ^ while all other modes are equal to zero 
identically. Then, we obtain from Eqs. (p6D p^: 

/i4 + 3 /i4 = —2.598 7 /ig , (30a) 
/is + 3 /is = 2.598 7 /il , (30b) 
(//, =0 forjV4,8). 

These equations describe the two-mode dynamical regime where both modes /i4(t) and 
/i8(t) (they form a conjugate pair) are connected with each other by the force interaction. 
Indeed, if the first oscillator, described by Eq. (|30|a) was excited (/i4(t) ^ 0),then the force 
2.598 7 /i4, acting on the second oscillator, appears in the right-hand side (rhs) of Eq. (pO|b) 
and, therefore, /is(^) cannot be equal to zero identically (/is(^) ^ 0). Conversely, the 
excitation of the mode figit) leads to the excitation of the mode /i4(t) because of the force 
term —2.598 7 /ig in Eq. (^a). Therefore, the dynamical regimes /i4(t) ^ 0,/is(t) = and 
/i4(t) = 0,fig{t) ^ cannot exist, and we may speak about the force interaction between 
modes f^iit) and fisit) (for more details, see 0, p!0|). 

The existence of the exact equations (|30| ) for the FPU-a chain with = 12 means 
that there exists a two-dimensional subspace 5*4 s of the twelve-dimensional configuration 
space spanned by the basis vectors (p^^ and (p^ (these vectors are associated with the modes 
/i4 and /is) which is dynamically invariant. Indeed, the trajectory corresponding to the 
dynamical regime (|30|) in the configuration space remains in the subspace S'43 for every 
time t. 

We do not discuss the problem of finding the analytical solution of Eqs. (|30|) in the 

present paper, and only want to note that these equations describe the two-dimensional 

^°Let us remember once more that the mode /ig is excluded from our consideration. 

^-'^ Hereafter, we keep in the coefficients of the dynamical equations only three figures after the decimal 

point. 
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bush with symmetry group G2[3a] as a dynamical object. 

Now let us consider the bush with G3[4a] from Eqs. (^). Assuming that /i3(t) ^ 0, 
IJ,G{t) ^ 0, /i9(t) ^ 0, and all other modes are equal to zero, we obtain from (|26D the 
following dynamical system with only three degrees of freedom: 

/i3 + 2/i3 = -47/i9/i6, (31a) 
fie + Aug = 27(/i3-/i9), (31b) 
/i9 + 2/i9 = A-ffisfiQ. (31c) 



Note that modes /i3 and fig are conjugate. They have the identical symmetry groups 
6*3 [4a], while the mode /ig possesses higher symmetry Gipa] D G3[4a]. Either mode fi^ 
or fig can be the root mode, and /is is the secondary mode. Indeed, if we excite the mode 
/i3(t) at the instant t = using the appropriate initial conditions, the force 27/i| acts on 
the mode /xg according to Eq. (pl|b) and, therefore, this mode turns out to be excited as 
well. (We can say this in another way: Eq. (|3T[b) cannot be satisfied in the opposite case 
fisit) ^ 0, fieit) = ). On the other hand, it follows from Eq. (|31|c) that figi^t) will be 
also excited when fisit) ^ and fj,e(t) ^ because of the force A'yfi^fiQ in the right-hand 
side of this equation. 

Thus, the excitation of the root mode ^^{t) brings about the excitation of the secondary 
mode /ig, but the opposite statement is not correct as a consequence of the specific asym- 
metric structure of bush dynamical equations with respect to the root and secondary 
modes. Indeed, the excitation of only the secondary mode /ig does not produce forces 
acting on modes /i3 and fig, because the forces i'jfiQfig and Ajfi^fis are equal to zero if 
Hsit) = and figit) = 0. Note that the previous geometrical analysis based on the con- 
sideration of the symmetry groups of modes (the symmetry of /ig is higher than that of 
/i3, fig) leads to the same conclusion which was obtained above by the studying of the 



structure of the dynamical equations (0). This is the general result proved by us in [|T^ 

(see Theorem 1): dynamical arguments must confirm the geometrical (group-theoretical) 

argumentsQ. The similar consideration of the bush with symmetry 6*4 [6a] leads to the 
Indeed, it follows as a consequence of this Theorem that when the group-theoretical method per- 
mits the existence of a nonzero invariant vector and, therefore, the existence of the secondary mode 
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following dynamical system with five degrees of freedom describing the evolution of modes 



/i2, /i45 /^65 /^8 and ^iq: 

/i2 + /^2 = 7 [-1-732 /i4/iio - 3.466 /ieAig] , (32a) 

/i4 + 3/i4 = 7 [ 0.866 /i^- 2.598 /i^- 3.466 yU6/xio], (32b) 

fid + 4/i6 = 3.4667 [ /i2/U4 - /U8/iio ] , (32c) 

/i8 + 3/i8 = 7 [ 2.598 /i^- 0.866 /i?o + 3.466 /is/ie], (32d) 

/Uio + Aiio = 7 [1.732 /i2/i8 + 3.466 /i4/i6] . (32e) 



Here conjugate modes /i2, /iio are root modes because they have the lowest symmetry 
G4[6a]. Conjugate modes /i4, fig with the group 6*2 [3a] D G4[6a] and mode /ie with 
Gi[2a] D 6*4 [6a] are secondary modes. This geometrical fact can be again checked by 
the appropriate consideration of the structure of equations (0). Indeed, the excitation 
of the mode ^2(t) implies the appearance of the mode /i4 because of the force 0.866 7/i| 
in Eq. (p2[b), and then the mode /i8 appears due to the force 2.598 7/i| in Eq. (|32|d). The 
presence of modes /i2, /i4, in the dynamical process leads to the excitation of the 
mode /ig, due to the force 3.466 7/i2/i4 in Eq. (^c), and to the excitation of the mode 
yUio because of the force 7(1.732 /i2 /i8 + 3.466 /i4/i6) in Eq. (|32|e). Thus, the dynamical 
equations (^) describe the evolution of the five-dimensional bush with symmetry 6*4 [6a]. 

The last bush from Eqs. (|20|) is the trivial one, because it contains all 12 modes and 
describes the dynamics of the FPU-a chain with twelve particles in the general case. 

Thus, we found all possible bushes of modes for the case = 12 generated by the 
subgroups of the parent group Go = ^ of pure translations, and it is clear that anyone 
can proceed in the similar way for considering the case of an arbitrary A^. 

In conclusion, let us consider the thermodynamical limit N 00 for the equation of 
motion (|27|) in the real modal space. 

corresponding to a given root mode, there also exist the invariants of the parent symmetry group which 
provide the force terms in the rhs of the dynamical equations leading to the excitation of this secondary 
mode. 
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2.4 Integro-difFerential equation for the FPU-a chain in the 
Umit N — > oo 

As usual in the physics of crystals, it is convenient to introduce the wave number k instead 
of the mode number j for the case ^ 1, 

271 

^=^J, j=0,l,2...,Ar-l. (33) 

All values of k belong to the interval [ 0, 2n) which corresponds to the primitive cell of 
the reciprocal lattice with the period b = where a is a period of the one-dimensional 
lattice associated with the FPU chain in equilibrium (we assume a = 1). The density of 
the fc-points from Eq. (|33D increases with increasing of A^. These points are distributed 
continuously and uniformly on the interval [ 0, 2tt) in the limit A^ — > oo. 

The number of terms on the right-hand side of each Eqs. ( ^71) increases as A^ — ^ oo, 
and we can transform the sum into a definite integral, assuming ^^(t) to be a smooth 
function of its index k. As a result of this procedure, the following integro-differential 
equation instead of the system of ordinary differential equations ( P7| ) can be obtainecQ: 

2tt 



z{x, t) + uj'^{x)z{x, t) = j z{y, t) \\^z{x + y,t) + z{-x - y, t)] K+{x, y) (34) 

z{x -y,t)- z{-x + y, t)] A'_ (x, ?/)| dy , 





+ 



X 

u;(x) = 2 sin(— ), A'-|-(x, y) = 2 sin(y) ± sin(a;), < x < 2tt . (35) 



where 



The function z{x,t) was introduced in Eq. ( |5^ instead of i^fc(t) and, therefore, z(x,t) 
signifies the value of the mode with wave number x at time t. Thus, we use different 
notations for the wave number in the discrete case {k) and in the continuous case (x). It 
is easy to verify the correctness of Eqs. (0, ^5]) by replacing the integral on the right-hand 
side of Eq. ( |3^ ) with its Darboux sum and using the following relation between z{x,t) 
and Vk(t)\ 



z{x,t) = —— Uk(f) , {x = k). (36) 

/TT 
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Note that the coefficient a can be removed from Eqs. (P7|) as well as from Eq. (|34|) by the appropriate 



scale transformation of the functions Vkit) and z{x,t)^ respectively. In Eq. (34) we also imply that the 
variables ±a; and zty are reduced to the interval [0, 27r) with the aid of cyclic conditions, x = x ± 27r and 
y = y±2-K. 
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Such a procedure returns us to the discrete model (P^). 

Note that Eqs. ( |5^ , ^) correspond to the continuum hmit for the FPU-a model in 
modal space, while the well-known Korteweg-deVries equation (KdV) corresponds to the 
continuum limit in ordinary space [ffl. Besides the form of above mentioned equations, 



there exists a principal difference between these two equations. Indeed, Eqs. (|3^ , |35|) , in 
contrast to the KdV-equation, is valid not only for long waves, but for short waves as 
well. We can make use of this fact for considering bushes of modes which often correspond 
namely to sufficiently short waves. The application of the integro-differential equation (|3^) 
for studying some problems of the FPU-a dynamics will be considered elsewhere. 



2.5 On the difference between bush structures for FPU-a and 
FPU-/? chains 

Let us compare our results, which are based on the group-theoretical method (Sec. 2.2), 
with those by Poggi and Ruffo , which are based on specific dynamical properties of the 
FPU-/3 chain. Naturally, taking into account concrete properties of a given mechanical 
system, we can obtain some additional selection rules as compared with general symmetry- 
related results. Poggi and Ruffo found for the FPU-/5 model five one-mode solutions 
corresponding to the following numbers of modes 

. ^ N N N 2N 3N , „^ 

J 4 ' 3 ' 2 ' 3 ' 4 ^ ' 

and two types of solutions associated with the pairs of modes {^, ^} and {y, ^}. Our 
results differ from these, and we want to discuss the problem using the example of the 
chain with = 12 considered in Sec. 2.2. 

According to Eqs. (^0|) there exists only one one- dimensional bush B[2a] associated 
with the excitation of the mode with j = ^ = 6 (remember that we, as well as Poggi 
and Ruffo, neglect the mode fio). But the excitation of the mode with j = y = 4 leads 
to the appearance of the two-dimensional bush B[3a] with j = 4,8 (modes /i4 and /ig 
are conjugate). Thus, in contrast to the FPU- (3 model, the one- dimensional dynamical 
regime with only fi^ (or fig) does not exist. Moreover, this conclusion is confirmed also by 
the (iynamzca/ consideration for the FPU-a model in the Sec. 2.3. The difference between 
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structures of the bush B[3a] for the FPU- (3 and for the FPU-a chains is brought about 
namely by their specific dynamical properties. 

Because of the same reason, the excitation of the mode fj^sU = x) leads, for the 
FPU-a chain and for the general well, to the excitation of the conjugate mode 

fJ-dU = ^) and to the secondary mode fiQ^j = y). We see, therefore, that some modes 
in the list (0), obtained by Poggi and Ruffo for the FPU-/3 chain, combine into the 
three-dimensional bush B[4a] for the FPU-a model, as well as for the arbitrary nonlinear 
chain. In other words, the three-dimensional symmetry-determined bush -B[4a], being 
an indivisible nonlinear dynamical object for FPU-a model, can be split into smaller 
invariant manifolds (i.e., bushes of modes in general sense) by imposing the appropriate 
initial conditions. Indeed, the one-dimensional bush with j = ^ can exist in the FPU-/5 
model according to Eq. (0), but it does not exist in the FPU-a model because it follows 
from Eq. ([31|b) that excitation of the mode ^^{j = ^) always leads to the excitation of 
the (secondary) mode fiQ^j = y). 

Let us stress once again that our results, obtained by the group-theoretical method, are 
general: they are correct for any nonlinear monatomic chain|^. Nevertheless, the specific 
character of interatomic interactions can reduce the dimensionality of bushes found by 
this method, as we have just seen comparing the structures of bushes of modes for FPU-a 
and FPU-/? chains. 

2.6 Bushes of normal modes associated with dihedral symmetry 
of nonlinear chains 

Introducing the concept of bushes of normal modes, we started with a certain symmetry 
group Go (parent group) of the considered mechanical system and basis vectors of its 
irreducible representations. Each given bush corresponds to some subgroup Gk of the 
group Go(Gfe C Go). 

In previous sections, the group T = E,a,a^, . . . ,a'^~^ (see Eq.(4)) was chosen as 
the symmetry group Go of the monatomic chain in equilibrium. But the full symmetry 
of such a system is higher, because space inversion i with respect to the center of the 
^''in some sense, they are also valid for multiatomic chains (see Sec. 2.6). 
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finite monatomic chain also is a symmetry element, as well as the product ■ i of the 
translation a''{k = 0, 1, . . . , — 1) with the inversion i. Therefore, the full symmetry 
group D (dihedral group) of the monatomic chain possesses twice as many elements as 
the translational group T which is its subgroup of index 2. We can write the group D as 
the direct sum of two cosets with respect to the subgroup T: 

D = T + T-i. (38) 

It is convenient for our purpose to rewrite this equation in the form 

Go = Go + Go- 1 (39) 

where Gq = T and Gq = D. 

We will discuss the connection between bushes of modes associated with the group 
Go and those associated with its supergroup Go- The irreducible representations of the 
groups Go and Go will be denoted by Tj and Tj, respectively. 

Up to this point we discussed bushes for monatomic chains in terms of the irreps Tj 
of the translational group Go = T, although the full symmetry group of these systems is 
Go = D. Why did we proceed in such a manner and how is the list of bushes of normal 
modes modified if we consider the complete symmetry group Go instead of its subgroup 
Go? We give the brief answer to this question below, while the detailed consideration will 
be given in another paper devoted to multiatomic chains. 

The group Go = -D is a non-Abelian group since some its elements do not commute 
with each other (for example, i ■ a = ■ i). As a consequence, the number of classes of 
conjugate elements of this group is less than the total number (2A^) of its elements and 
some irreps Tj are not one-dimensional. 

Let us discuss the monatomic chain with A^ = 12 particles in more detail. The 
appropriate group Go = -D has nine classes of conjugate elements and, therefore, nine 
irreducible representations Tj. Four of these irreps (Fi, Fs, r4) are one-dimensional 
and five other irreps (F5, Fg, F7, Fg, Fg) are two-dimensional. 

To determine a given irrep F^ of the group Go defined by Eq. (|39|) , we must write 
matrices of this representation for elements of the first coset Gq and of the second coset 
Go • i- They form the first half vf^ and the second half f ^-^^ of the irrep f j, respectively.^ 

Obviously, v'^^ is the restriction of the irrep T j of the group Gq onto the subgroup Gq . 
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fj.^^ for one-dimensional irreps coincides with either irrep Fq or irrep Tq of the group Go 
(see Eq. ( |¥TD below). T^^^ for two-dimensional irreps can be written as a direct sum of two 

conjugate irreps of group Gq, for example, = 1 . All matrices of the second 



half rf of the irrep Tj are products of the appropriate matrices of F^^"* with the matrix 

corresponding to the inversion (i). For one-dimensional irreps f ^ this matrix is equal to 

/O 1\ 

(-1-1) or (-1), while for two-dimensional irreps F,- it is equal to . Therefore, the 

\l o) 

irrep F5 has the following form^ 

rii)^( (1 
This view of the irrep Tj must be compared with the definition (^) of group Go', it allows 
one to determine the explicit form of all matrices of the both cosets from Eq. (pOf). 



(40) 



Using the above notation we can write all irreps Tj {j = 1, ...,9) of the group Go as 
follows: 

fi = Fo + Fo-(l); f2 = Fo + Fo-(-l); 

f3 = F6 + F6-(l); f4 = F6 + F6-(-l); (41) 

■Ffc_4 \ /Ffc_4 \ / 1 

F,= I + I , A; = 5,6,7,8,9. 



Note that irreps Ffc_4 and Fig-fc in Eq. (PD are mutually conjugate for any k: 

Ti-Fn, F2-F10, F3-F9, F4-F8, F5-F7. (42) 

As was already discussed, bushes of normal modes correspond to different subgroups 
of the parent symmetry group of the considered mechanical system in equilibrium. For 
monatomic chains, bushes associated with the parent group Go were considered in Sec. 
2.2. Now we want to discuss bushes associated with the parent group Gq. 

There are 34 subgroups of the group Go for the chain with = 12 particles. Note 
that since Go C Go, all subgroups of Gq are also subgroups of Go. Therefore, all bushes 
found in Sec. 2.2 are valid not only for the parent group Gq, but for the parent group 
Gq, too. In addition to these bushes, some new bushes of normal modes can appear when 
we consider Go as the parent group. To find them we must examine those subgroups of 



Henceforth, we will omit zero elements of all matrices. 
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Go which contain some elements of the second coset (Go • i) from Eq. (^). Obviously, 
all these new subgroups can be constructed by adding certain elements of Go ■ i to the 
subgroups of the group Gq (in particular, to the trivial subgroup {E}). Let us consider a 
concrete example. 

The subgroup {a^,i} of the group Go is determined by two generators and i (here- 
after we write down in braces only generators of the appropriate group instead of all its 
elements). This subgroup is obtained from the subgroup {d^} of the group Go = T by 
adding inversion i. As was shown in Sec. 2.2, the three-dimensional vibrational bush 

5{a^} = /i6(t)¥'6 + l^3it)(fi3 + lX9{t)<p9 (43) 

corresponds to the subgroup {d^} where three functions fie(t), fi3{t) and fig(t) are inde- 
pendent of each other. Thus, three one-dimensional irreps Fg, and Fg contribute to the 
bush B{d'^} associated with the parent group Gq = T (namely their basis vectors are 
contained in Eq. (^3]). 

In terms of irreps of the parent group Go, the same bush is formed by basis vectors 
of one-dimensional irrep F3 and two-dimensional irrep F7 (note that restrictions of these 
irreps to the subgroup {d^} give us the above mentioned irreps Fg, F3 and Fg of the group 
Gq = T). The additional generator i which appears when we pass from the subgroup 
{d^} to the subgroup {a^, i} introduces the matrix (1) for the irrep F3 and the matrix 
for the irrep F7, and the new bush B{d,i} must be invariant with respect to 
these matrices. This invariance condition leads to the following transformation under the 
action of i: 

IJ'&it) I^Git), Ai3(^) l^git), yUg(t) /i3(t) . (44) 

As a result, we obtain the restriction yU3(t) = /ig(t) on previously independent functions 
/i3(t) and fig{t). In other words, the old three-dimensional bush B{d'^} from Eq. ( |45| ) 
transforms into the new two-dimensional bush 

B{d\ 1} = /i6(t)c^6 + /^3(t) [V3 + ^9] (45) 

Let us emphasize that the transition from the parent group Go = T to the parent 
group Go = D results in new superpositions of one and the same set of normal modes 
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'^j (j = 0, 1, ... , 11) when we find the appropriate bushes of modes and, therefore, we 
can compare bushes associated with different parent groups. 

Thus, all bushes associated with the group Gq are also bushes associated with the 
group Go, but there appears a set of new bushes for Go because of additional restrictions 
on time- dependent functions for some old bushes. The detailed analysis of this 

phenomenon will be published elsewhere, while we give the appropriate results in Table 2 
of the present paper. 

The dimensions of the bushes are indicated in the first column of this table. The list 
of the appropriate irreps Tj in square brackets in column 2 corresponds to a certain bush. 
Such a list determines the bush uniquely, since a quite definite basis vector Lp^ corresponds 
to the irrep Vj of the group Go. 

For bushes associated with the parent group Gq = D (column 3) we give in square 
brackets not only individual irreps Tj (of the group Gq = T!) but also the pairs Fj — Tj of 
conjugate irreps (the irreps of such pairs are connected by hyphens). Only one function 
fii{t) corresponds to each pair Fj — Tj. The above discussed bushes (|43|) and ( ^Sf ) can be 
written in this notation as follows: 



The relation between fisit) and ^dit) for the bush B{a^, i} (see Eq. (|^)) is very simple 
(yU3(t) = fig{t)), while for some other bushes the similar relations can be more complicated. 
Namely because of this reason we do not define concretely the relation between functions 
fii{t) and fJ'jit), and write simply Fj — Tj, indicating by the same token only that such 
a relation does exist and, therefore, instead of two independent functions fii{t) and /ij(t) 
only one function /ii(t) corresponds to the pair Fj — Fj. 

Moreover, each square brackets in the third column of Table 3 can correspond to 
several bushes which differ from each other by different relations between fii{t) and /Uj(t) 
associated with pairs Fj — Fj. By the way, this is one of the causes why the number of 
bushes in Table 3 is less than the total number (34) of subgroups of the group Gq = D. 
There exist also another cause of the above phenomenon. Indeed, there are no vibrational 
bushes for some subgroups of the parent group Gq = D. As examples we can point to 



s{an = [r6,r3,F4], 



(46) 



(47) 
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subgroups {a,i} and {a?,di}. Finally, it is possible that adding a new generator does not 
change the appropriate old bush (for example, B{a^,i} = B{a'^}). 

To conclude the discussion of relations between bushes associated with parent groups 
Go and Go, which are indicated in Table 2, let us say that besides five bushes associated 
with Go there exist seven new types of bushes in accordance with their classification 
under the parent group Gq: 2 one-dimensional, 2 two-dimensional, 1 three-dimensional, 
1 five-dimensional and 1 six-dimensional. 

Let us stress that for some multiatomic chains inversion is not a symmetry elementp^ 
and Go = T is the full symmetry group of these mechanical systems. In such cases we only 
need to use Go = T as the parent group for finding normal modes. For monatomic chains 
and for multiatomic chains with inversion, using Go as a parent group provides more rich 
information about possible bushes and, therefore, is preferable. In most sections of the 
present paper, we restrict ourselves to studying only bushes associated with the group 
Go = T, since this case is more general and more simple. 

It is very essential that our list of bushes given in Sec. 2.2 (see also the second column 
of Table 2) is valid for every one-dimensional structure with the symmetry group Go = T, 
independently of how many particles are in its primitive cell: the difference appears only 
in the form of basis vectors^ but not in the bush structure. Let us discuss this point using 
as an example the bush with subgroup G3[4a] from Eq. (^) which was denoted in the 
present section as 5{d^} (Eq. (^3])). 

The bush of such type can be excited in various chains, but the basis vectors c^^ 



from (|43| ) will be different for different structures of the primitive cell. In particular, the 

dimension of the vectors c^^- and X{t) is equal to A; ■ if there are k particles in the 

primitive cell. Moreover, the dimension of the bush also can change, because, in general, 

basis vectors ^p^ depend on a number of arbitrary constants and this number is equal to 

the number of copies of the appropriate irreps in the decomposition of the full vibrational 

representation of our mechanical system into its irreducible parts [10]. 

Thus, the group-theoretical methods give us a possibility to divide the problem of 

finding all types of symmetry-determined nonlinear vibrations into two independent parts: 
^^In terms of crystallography, we can say that it is possible for a given crystal to not have inversion 
even though its Bravais lattice does. 
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1) obtaining the list of bushes using group-subgroup relations only; 

2) taking into account the concrete structure of the considered mechanical systems. 
We will discuss bushes of normal modes for multiatomic chains in more detail in a 

special paper. 

3 Stability of bushes of normal modes 

Some aspects of the stability of exact solutions for the FPU-/? chain were discussed by 
Poggi and Ruffo in . They used the standard method based on the idea of the lineariza- 
tion of appropriate dynamical equations in the vicinity of a given solution. We also make 
use of this method, as well as of a direct numerical study of the dynamical equations for 
the FPU-a chain in the ordinary and in the modal space. 

The loss of the stability of bushes for small and for large N can occur in different 
ways. It is convenient to discuss these cases separately. 

3.1 Stability of bushes of modes for small 

Let us consider the problem of bush stability using as an example the FPU-a chain with 
N = 12. All numerical data on thresholds of the bush stability, given in the present 
section, correspond namely to this case. 

Firstly, we discuss the stability of the one- dimensional bush 5 [2a] with the symmetry 
group G = [2a]. It contains only one mode i^Qit). Therefore, we must consider all other 
modes (j ^ 6) in Eqs. (pTf ) to be equal to zero in the exact solution corresponding to the 
given bush. Bearing in mind the intention to study the bush stability, we suppose that 
these "sleeping" modes are not exactly equal to zero, but are sufficiently small quantities 
(for example, at about 10^^° for the numeric investigation). Then the terms of Eqs. (|27| ) 
differ from each other by the order of their smallness: 

I ^j^6 I > I ^i^k I for i,j,ky^6. (48) 

As a result, we can keep only those terms in the right-hand sides of Eqs. (^Tf ), which 
contain the mode z/6(t). Thus, these equations can be reduced for A^ = 12 to the following 
form (hereafter all numerical coefficients are given up to the third figure after the decimal 
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point): 



i)i + ufui = — 1.155 a Z/5Z/6 , (49a) 

i>2 + ujIu2 = —2 a U4^uq , (49b) 

i>3 + cuius = —2.309 a u^uq , (49c) 

i>4 + ujlh'4 = —2 a z/2t'6 5 (49d) 

i>5 + cUgZ/s = — 1.155 a z/iz/g , (49e) 

h + ujlu6 = 0, (49f) 

i>7 + cUyZ/y = 1.155 a z/iiz/g , (49g) 

i)s + ujlus = 2 a uioUq , (49h) 

i)g + cUgZ/g = 2.309 a u^uq , (49i) 

i>io + cj^gz/io = 2aust^e, (49j) 

1)11 + ujI^uu = 1.155 a ujUq , (49k) 

^2 = 4sin2f^') , j = l,2,...,ll. (50) 



where 

.127 

The system ( ^9]) of eleven equations sphts into seven independent subsystems. Moreover, 



it follows from the equation (HSlf)^^ 

^6 + ^^6^6 = {uji = A) (51) 
that the root mode z/6(t) vibrates in harmonic manner , 

z/6(t) = Acos(2t + (5), (52) 
and, therefore, we can substitute u^it) in this form into all of the other equations in (^9]). 



Thus, we obtain 



i)i + cului = — 1.155(ay4)z/5 cos(2r) , (53a) 
i>5 + cUgZ/s = — 1.155(ay4)z/i cos(2r) , 
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Let us emphasize that this is the exact equation of the one-dimensional bush B[2a] (see also Eq. ( P9| ) ) . 
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i)2 + uJ2^2 = — 2(aA)z/4 cos(2r) , 
1^4 + C1J4Z/4 = —2{aA)u2 cos(2r) , 



(53b) 



u^ + uj^us = -2.309(aA)z/3Cos(2r) , (53c) 

i>7 + co'yZ/7 = 1.155(aA)z/ii cos(2r) , (53d) 

+ cul^uii = 1.155(a74)i^7 cos(2r) , 

i>8 + tUgZ/g = 2(ay4)z/io cos(2r) , (53e) 

i^w + ^^IqI^w = 2{aA)us cos(2r) , 



ug + ulug = 2.309(aA)z/9Cos(2r) . (53f) 

The initial phase 6 was removed from the above equations by introducing the new time 
variable 

T = t+^-. (54) 

Individual subsystems (^) describe different ways that the stability of the original bush 
B[2a] can be lost. Indeed, we obtain conditions of the loss of stability with respect to the 
modes i^s^t) and I'g^t) from Eqs (|53|c) and (p3|f), respectively, which can be easily reduced 



to the standard form of the Mathieu equation ||18|| : 

i/+[a-2gcos(2r)]y = 0. (55) 

The set of stable and unstable regions correspond to this equation in the plane (a, q) of 
its pertinent parameters. Using the stability chart of the Mathieu equation in the same 
manner as in |^ we obtain, for the considered way of the loss of stability, that the bush 
i?[2a], describing by the mode z/6(t), is stable for 

\aA\ < {aA)^ = lM9. (56) 

Moreover, it can be stable for larger values of | a A |, as well (for example, for the 
interval 16.465 <| a A \< 16.474), but in this paper we study only the first (basic) zone 
of stability for all bushes of modes. 
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Appearance of the modes and z/g means that the original one- dimensional bush 
5 [2a] enlarges up to the three-dimensional bush B[4a] which embraces the root modes 
z/3,z/9 and the secondary mode z/g [see Eq. (pO|)]. Let us rewrite Eqs. (^9|c), (^9|f) and (^), 



describing the dynamics of this bush: 

i)^ + 2z/3 = -2.309 a v^^uq , (57a) 
i>6 + 4z/6 = 0, (57b) 
uq + 2z/9 = 2.309 a ugUQ . (57c) 

Eqs. (^) represent the linearized equations (with respect to the "sleeping" modes and 
uq) of the exact equations of the bush i?[4a] in terms of real modes^ 

2 8a 

1^3 + ^3^3 = 7f^^^3i^6, (58a) 

V iz 



2 4a 



i)e + ujQi^e = --7f^ i^3 -^9)^ (58b) 



2 , ^ , 

UQ + ujgUg = -y=UQUe. (58c) 



In contrast to the equations (|3T|) for the bush _B[4a] in terms of complex modes /is, /Xg, /^g, 
Eqs. (|58|a) and (|58|c) are independent from each other. Nevertheless, the parametric 
excitation of the mode i^slt), brought about by the mode z/6(t) = Acos{2t + 6), leads 
simultaneously to the excitation of the mode lygit). Indeed, the same condition of the 
loss of the stability corresponds to both equations (p7|a) and (p7|c), because they can be 
converted to the identical form using the transformation (|5^) with the initial phase shifted 
by TT. 

Equations (§^b) for modes z/2(t), z^4(t) and (^e) for conjugate modes z/io(t), J^sit) 
describe the transition from the original bush B[2a] to the embracing bush 5 [6a] because 
of the parametric resonance with the mode I'eit) [see (0)]. Note that Eqs. (|53|b) and (|53|e) 
can be converted to the same form using the above mentioned transformation reducing 
Eqs. (p7|c) to the form of Eqs (|57|a) and, therefore, we may discuss equations for modes 
z/2(i) and 1^4 (t) only. 

It is interesting to stress that we cannot consider the parametric excitation of the bush 
-B[3a] by the mode z/6(t) belonging to the bush i?[2a]. Indeed, according to Eq.(^) the 



29 



Compare these equations with those in terms of complex modes (^l|). 
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two-dimensional vibrational bush 5 [3a] consist of two modes, z/4 and ug, whose symmetry 
groupQ is G = 3a. On the other hand, the new bush B, appearing as a result of the 
loss of stability of the old bush i?[2a], must include the last bush and, therefore, must 
contain the mode uq. Thus, the above mentioned loss of stability can induce only bushes 
5 [6a] and i?[12a], but not the bush B[3a] (see Eq. (^OD). This geometrical fact can be 
confirmed by the following dynamical arguments. 

We find from Eqs. (|53| ) that the active mode ueit) (see Eq. (|52D) can excite the following 
pairs of the sleeping modes ui, z/5 (Eq. (|53|a)), z/2, 1^4 (Eq. (p3|b)), z/y, uu (Eq. (|53|d)) 
and Us, ^10 (Eq- (His))- AH modes in the pairs (z/i, z/5) and (z/7, z/n) possess the same 
symmetry group G = 12, while the modes from each other pair, i.e. (1^2, 1^4) and (z/g, z^io), 
have different symmetry groups: the group G = 3a corresponds to z/4, z/g and G = 6a 
corresponds to z/2, z/iq. Since both modes of each of the above mentioned pairs must 
appear simultaneously (this can be seen from the structure of Eqs. (|53|a), (|53|b), (^3]d) 
and (p3^), respectively), the bush 5 [3a], in its pure form, cannot be excited in connection 
with the loss of stability of the bush i?[2a]. 

Indeed, as was just stated, according to, for example, equation (|53| b), the excitation 
of the mode z/4 with G = 3a leads necessarily to excitation of the root mode z/2 {G = 6a) 
of the larger hush i?[6a]. Thus, equations (p3|b) 

i)2 + ujliy2 = -2(aA)z/4Cos(2r) , (59) 
1/4 + (^4 z/4 = — 2(aA)z/2 cos(2r) 

describe the transition from the original bush i?[2a] to the bush 5 [6a] [root modes z/2(t), 
z/io(t)], when {a A) reaches the appropriate critical value {aA)^. Direct numerical compu- 
tations for Eqs. (|5^) lead to the following result: 

(aA)c = 1.049. (60) 

Comparing the critical values (|56|) and (|60|) of (aA) providing the transition to the bushes 
B[4a] and -B[6a] respectively, we see that these bushes are excited by the mode z/6(t) at 
the same threshold (up to the numerical accuracy). 
30 See Table 1. 
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Similar to the above analysis, we can consider the stability of the bush B\2a\ with 
respect to interactions with the root modes of the bush i?[12a]. There are four such modes: 
z/i, 1/5 and z/7, Vii conjugate to them. Nevertheless, it is sufficient to take into account only 



two equations (|53|a) for Pi and z/5 because Eqs. (|53|d) for z/7, z/n can be transformed to the 
form (|53| a). The critical value {aA)c which corresponds to the transition 5 [2a] — > 5 [12a] 
also coincides with the above values (|56|) and (|60|) for transitions 5 [2a] i?[4a] and 
B[2a] B[Qa]: 

(aA)c= 1.049. (61) 

We will return to this surprising coincidence in the next section of the present paper. 

Comparing all three variants [see Eqs. (0), (160|), (|6lD] of the loss of the stability of 
the original bush 5 [2a], which brought about its extension up to the bushes -B[4a], -B[6a] 
and -B[12a], respectively, we find that as the amplitude A of the mode ^^{t) is increased, 
the transition i?[2a] i?[12a] must occur because the bush 5 [12a] embraces both bushes 
i?[4a] and i?[6a]. This result, obtained from approximate Eqs. (p5D, was confirmed by 
direct numerical computation of the exact dynamical equations for the FPU-a model. 
We wrote two variants of the appropriate computer programs: the program PI solves 
Eqs. ( p3D in ordinary space and then decomposes the vector X(t), defined by Eq. @, into 
the set of modes, while the program P2 solves Eqs. (pT]) for the FPU-a chain in modal 
space. Both programs give the same estimation for the value of {aA)c: 

1.04 < < 1.05 , 

and confirm the above conclusion that the bush i?[2a] really transmutes into the bush 
i?[12a] as a result of the loss of its stability. 

Below we give the critical values of [aA) for the loss of stability of all other bushes 
for the FPU-a chain with N = 12: 

5[4a] ^ S[12a] : ^ 1.05, (62) 

B[3a] B[6a\: {aA)c^ 0.60, (63) 

B[6a] 5[12a] : (aA)c ^ 1.29. (64) 

Let us pay attention to the loss of the stability of the bush i?[3a]. As follows from 
Eq. (|63D , this bush transmutes at {aA)c = 0.60 into the embracing bush i?[6a]. This 
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transformation is accompanied by the lowering of symmetry (G[3a] D GfGa]) twice. The 
new bush S[6a], arising as a result of this process, exists as an individual stable object 
for the following interval of its own root mode 

0.60 < av4 < 1.29. (65) 

In turn, the bush B[Qa\ loses stability for [aA) > {aA)^ = 1.29 and transmutes into 

the largest bush B[12a] containing all modes uj [j = 1, 2, . . . , 11), also with lowering of 

symmetry twice {G[6a] D G[12a]).Q 

Thus, we have considered the stability of all possible bushes for the FPU-a chain with 
= 12 particles.0 All other cases with sufficiently small can be studied by similar 

methods. But some new phenomena appear when A^ becomes very large, and we consider 

them in the next section. 

3.2 Stability of bushes of modes for large N 

Let us consider the stability of the bush B[2a] for an arbitrary even value of A^ (its stability 
for the case A^ = 12 was discussed in Section 3.1). This bush consists of only one mode 
with j = N/2 and the appropriate dynamical equation 

+ ^N/2^N/2 = (66) 

can be obtained from Eqs. (|27|) assuming that all modes different from i'N/2{t) are equal 
to zero. This is equation of the harmonic oscillator with ujn/2 = 2 [see Eq. (|25|)], and we 
can write its solution in the form: 

z/Ar/2(r) = Acos{2t) . (67) 
■^^Transition from one bush Bi to another bush B2 with higher dimensionahty (with more degrees 
of freedom) as a result of the loss of stability of Bi is always associated with spontaneous breaking of 
symmetry of the original dynamic regime described by the bush Bi . This is one of the central conclusions 
of the bush theory since every bush possesses its own symmetry group In connection with this, let 

us note that some aspects of the loss of permutational symmetry of the vibrating chains are considered 
in (§]. 

In this paper, we consider stability of the bushes associated with the translational parent group 
Go = T only 
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Linearizing the system ( JI7\ ) near the exact solution (^) leads to the following approximate 
equations: 

Uj + ujjUj = -— = sm j , (68) 



^\ . _ . . ^ _ 1 ^ 



= 4 sin2 ( ) , J = 1, 2, . . . , - - 1, - + 1, . . . , iV - 1 . 



Note, that right-hand sides of these equations contain only the term (the 
coefficients before similar terms Z/7V/2 ^N/2+j turns out to be zero). It is easy to see that 
the system (|68D splits into a number of subsystems containing one or two equations only. 
Indeed, the mode Uj in (|68D is connected with the mode u-, where j = y — j, and vice 



versa, the mode is connected with i^n_j = uj. Therefore, we have for j = 1, 2, . . . , y — 1 
in Eq.(^) the following pairs of equations which are independent {toiu. all other equations: 

i)j + 4 sin^ (^^^ Uj = -7 sin (^^^ un_j cos(2r) , (69) 
i)N^^ + 4 cos^ (^^) UN_^ = -7 sin (^^) cos(2r) , 



7^^, .-l,2,...,(f-l). (70) 



where 



Here we made use of the explicit form of i'N/2(t) from Eq. (l67|). For j = y + 1, . . . , — 1 
we obtain from (^) the equations similar to (^), but with the changing of the sign of 
the constant 7. As was discussed in the previous section, this case can be reduced to 
Eqs. (^) by replacing r ^ r + |. (Note, that for j > ^5 the index j = y — j becomes 
negative and then we can use the cyclic condition j + N = j for obtaining positive values 
of j). When j = y — j, both modes and i^N_j reduce to the same mode un and (|69|) 
reduces to the pair of identical equations of the form 

i>jv + 2i/]v = —7 z/jv cos(2r) . (71) 

44 4 

Obviously, this equation can be transformed to the standard form of the Mathieu equation 
(^) with ?/(r) = i'n^t), a = 2,q = —7/2. 

Thus, the system ( |68|) of — 2 differential equations splits into the subsystems (|69|) 
or (^) according to the diagram shown in Fig. 1, where indices j corresponding to pairs 
of equations (|69|) are connected with the arches of different size. 
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It is convenient to rewrite Eqs. ( pD| ) as follows:^: 

X + 4 sin^(A;)x = 7 sin(2A;)y cos(2r) , (72) 
y + 4 cos^{k)y = 7 sm{2k)x cos(2r) . 

where k = ttJ/N, x{t) = Vjir), y^r) = i'n_j{t). Thus, studying the loss of stabihty of the 
bush i?[2a], brought about by its interactions with the modes z/Ar/4 and I'sN/'i is reduced 
to analyzing the Mathieu equation (|7TD , and the loss of its stability with respect to all 
other modes reduces to analyzing Eqs. (^). 

We begin to study the stability of the bush B[2a] by considering its interaction with 
the mode un/a [let us remember that B[2a] consists of only one mode i'n/2 described by 
Eq. (pTDJ. It is clear that B[2a] must be a stable object for very small values of a A and 
we are interested in the lower boundary Ac of the loss of its stability, because of the 
parametric excitation of the mode i^n/a, when the amplitude A of the mode increases 



from zero. This problem can be solved easily with the aid of the Mathieu equation (71). 
Using its well known stability chart and analytical formulas for boundaries of unstable 
regions in the {a,q) plane [0, and independently by direct numerical investigation of 
Eq. (^), we found the following threshold value 7c for the coefficient 7 entering into this 
equation: 

7c = 2.42332 . (73) 

Then from Eq. (|70D we obtain: 

aAc = 0.30292ViV . (74) 

Thus, the bush B[2a] is stable with respect to the interaction with the mode I'N/'i for 
I a A \< a Ac given by Eq. (|7^) . It is obvious that the same value of a Ac corresponds to 
the loss of stability of B[2a] because of parametric excitation of the mode ly^N/i- 

Now we must take into account interactions of the bush B[2a\ with all other modes 
z/j(t). As was already stated, studying the stability of 5 [2a] leads, in this case, to the 
equations ([72|). They are linear differential equations with periodic coefficients and, there- 
fore, the Floquet theory can be applied to them. Unfortunately, we did not obtain any 
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As was already noted, the sign before 7 in Eqs. (p3) is not essential. 
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exact solutions of Eqs. (|72p and can now give only results of the appropriate numerical 



analysis. Our computation of the multiplicators for the system ( [T^ ) as eigenvalues of the 
monodromic matrix reveals a surprising fact! Indeed, it seems that the critical value 7c of 
the constant 7 from (|69|) [see also Eqs. (0)], corresponding to the loss of stability of the 
bush i?[2a], must depend on the mode number j, since the coefficients of these equations 
depend explicitly on j/N. However, this is not true. We found that 7c does not depend 
on j/N, at least up to 10^^, and coincides with that of the Mathieu equation (|7l|): 

7c = 2.42332. 

This nontrivial fact means that the original bush B[2a] loses its stability with respect 
to all modes uj {j 7^ N/2) simultaneously, i.e., for the same value of the amplitude 
of the mode i'N/2{t)- In other words, all modes of an A^-particle FPU-a chain are ex- 
cited parametrically because of interaction with B[2a] when a A reaches its critical value 
aAc = 0.30292-\/iV and, therefore, B[2a] transforms at once into the bush B[Na] of 
trivial symmetry.^ 

The above results obtained by studying Eqs. (^,^) were verified with the aid of the 
program PI for solving the dynamical equations (^) for the atomic displacements Xi{t), 
i = 1,2, . . . , N. Indeed, for = 50,100,200, we obtained in such a way the relation 
aAc = 0.303-\/iV which is in the excellent agreement with Eq. (fril). Moreover, it was 



also found that all modes i'j{t),j = 1,2, . . . , N — 1 are present in the decomposition (jlj) 
of the configurational vector X(t) from Eq. (j^). This fact confirms the above discussed 
conclusion about the appearance of the bush B[Na] from the original bush B[2a\ when 
the amplitude A exceeds the critical value Ac from Eq. (|74D. 

Note that using the program PI for studying the stability of the bush -B[2a], we try 
N only up to 1024, but analyzing this problem with the aid of Eqs. (|69|) we try very 
great values of A^, such as = 10^, 10^ allowing j/N to approach sufficiently close to the 
"dangerous" zero value. 

Now let us consider the configuration vector X(i(:) from Eq. ( |14|) which determines 
displacements of all particles from their equilibrium states expressed in terms of the 
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We already met this phenomenon [see Eqs.(p6|),(|60|) and (|6l|)] for the special case N = 12. The value 
aAc — 1.049 from the above mentioned equations is obtained from our present formula aAc = 0. 30292 
for N = 12. 
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modes i^j{t). The bush B[2a] consists of only one mode z/jv/2 and, therefore, 

X(t)|^p^j = z/^/2(t)V'iv/2- (75) 

The components {'4'N/2)n of the N-dimensional vector i/'Ar/2 can be obtained from Eq. (|TB|): 

(V'w/2)n = icos(7m), n = 1,2, ...,Ar. (76) 
V iV 

Then we have 

1 

and 



V'jv/2 = c, where c = (-1, 1, -1,1,..., -1, 1) (77) 



Here, for simphcity, we choose such an initial condition for the excitation of the bush i?[2a] 
that 5 from Eq. (|5^ ) is equal to zero. It is clear from (^), (^) that the displacement 
pattern of the FPU chain, corresponding to arbitrary instant t, possesses the translational 
symmetry 2a (see Fig. 2), since displacements of every two particles 2a apart from each 
other are identical. 

The dynamical regime ([78| ) loses its stability when the amplitude A becomes larger 
then Ac determined by Eq. ([741): 

« ^ ViV. (79) 

a 

Substituting this value into Eq. (^) we obtain the cntica/ form X(t) for the considered 
bush B[2a]: 

X.(t) = ^ c cos(2t). (80) 

B\2a] a 

Note that Xc(t) from (|80|) does not depend on the number of particles (A^) in the FPU 
chain! It follows from (^) that amplitudes of all atomic displacements for critical Xc()f:) 
must be larger for lesser values of a and, formally, they tend to infinity when a — 0. 
Such a behavior reflects the fact that there are no interactions between normal modes in 
the limit a = and, therefore, the normal mode VN/2{i) becomes stable in this case for 
any arbitrary amplitude A. To avoid a possible misunderstanding, let us note that for a 
FPU-a system, considered as an abstract exact model, any distance a between particles 
in equilibrium is permitted and, as a result, atomic displacements can be large. From 
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this point of view, the presence a in the denominator in Eq. (^) must not be striking. 
However, when we consider the FPU-a model as an approximation to the real physical 
system, only sufficiently small atomic displacements are assumed because we must neglect, 
in such a case, the higher terms in the Taylor series (|lD for F(Aa;). We must also take into 
account that the value of the amplitude A of the oscillatory regime is restricted by the 
presence of a potential barrier which appears because the nonlinear term with a in the 
FPU-a Hamiltonian (^) can be negative with absolute value larger than the quadratic 
term (the height of this barrier and its distance from zero are proportional to l/a^ and 
I /a, respectively). 

There was a discussion on the loss of stability of the zone-boundary mode (ZBM) in 
pO| , |19|.Q In our terms, this problem is identical to the problem of the loss of stability 
of the bush -B[2a], since this bush consists of only one above mentioned mode.Q In the 
present paper, we do not want to give detailed comparisons of the results by different 
authors and restrict ourselves by the following short comments. 

In the paper |Jl9[ by Sandusky and Page, the interrelation between the loss of stability 
of the zone-boundary mode (ZBM) and existence of intrinsic localized modes (ILM's)P^ 
is discussed for different nonlinear monatomic chains. It was found that there can exist 
two different types of the loss of stability of the ZBM in the (/c2, k^, k^) chains.]^ One of 
these types is connected with the appearance of the ILM's, while the other, the so called 
"period-doubling" type[^, is not related to ILM's. 

Only the "time-doubling" type of the loss of stability of ZBM occurs when the coef- 
ficient k^ is sufficiently large, in particular, when k^ ^ 0, k^j^ = (this case corresponds 
to the FPU-a model). Namely this phenomenon is discussed in the present paper for the 
FPU-a chain, and we can compare the values of the amplitude threshold for the loss of 
stability of the ZBM (the bush _B[2a], in our terms) obtained by us and by Sandusky and 

■^^ We are very grateful for the referee of our paper for these references which we did not know previously. 
■^^ In some papers, the zone-boundary mode is called by the term "7r-mode" . 
•'^ The term "breathers" is usually used for ILM's at the present time. 

^^Here, fc2,fc3,/c4 are coefficients in front of the harmonic term (^2) and in front of the anharmonic 

terms of the third {k^) and fourth (^4) orders. 

This term means that the new mode which appears because of the loss of stability of the ZBM 

possesses a time period two times larger than that of ZBM. 
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Page. The following estimate was found in [19]: 

0.302 <| A3 |< 0.303, (81) 

where A3 = k?jA/k2 with A being the non-normalized amplitude of the ZBM. 
Taking into account Eg. ([2^) , we can express A3 in terms of our paper: 



A3 = aA/^/N, (82) 

where A is the mode amplitude with respect to normalized basis vectors (p!3D. Then, with 
the aid of the Eq. @ 

aAr = 0.30292ViV 



we obtain 

A3 = 0.30292 (83) 

This value of the threshold of the loss of stability of the ZBM is in the excellent agreement 
with estimation (81) obtained in the paper [19] for the case = 100. Our result (|83|) 



is valid for arbitrary value of and, moreover, as was already discussed in the text of 
the present paper, we found that all modes appear simultaneousl]^ when A3 exceeds the 
remarkable threshold (|83|) . This threshold was obtained by another method, as compared 
to [19], and with the higher precision. 

Thus, we confirmed the validity of the results on the loss of the ZBM stability for the 
FPU-q; chain by Sandusky and Page from the paper [19] in contrast to those obtained 
in [20,21] (the discussion on the causes of the incorrectness of the appropriate results 
reported in [20,21] can be found in [19]). 

In many papers, the energy density e = E{0)/N of the initial excitation, i.e. the 
energy E{0) per one particle, is considered as a relevant control parameter characterizing 
the loss of stability of different modes in nonlinear chains. Substituting Xc{0) from (|80|) 



into the Hamiltonian (p2D leads to vanishing the sum of anharmonic (cubic) terms, and 
we obtain the following expression for the threshold ec of the loss of stability of the bush 
B[2a]: 

e = 0.18325/al 

"^^We suspect that this phenomenon leads, as a final result, to equipartition of the energy between 
difierent modes of the FPU-a chain, but we did not study this fact. 

43 



Let us emphasize that this value does not depend on the number of the particles in 
the FPU-a chain. 

We want to note that unlike one-dimensional bushes, the energy density e is not a 
relevant parameter for characterizing bush stability in the general case. Indeed, it can be 
shown that if a given bush contains several modes, the threshold of the loss of stability 
depends not only on e but on the distribution of the energy among its modes at the initial 
instant. Because of this reason, we prefer to fix this initial energy distribution in a certain 
way, namely, we assume that the energy is completely localized only in the root mode of 
the bush and then we look for the threshold value Ac of this mode. 

In conclusion, let us note that the rigorous stability analysis of all bushes of vibra- 
tional modes cannot be fulfilled on the basis of the Floquet theory^, because, in general, 
dynamical regimes corresponding to them are not periodical. Indeed, a given bush can 
contain modes with incommensurable frequencies (for example, the bush i?[4a] contains 
the modes with the harmonic frequencies tus^g = \/2 and ujq = 2). 

Now let us compare our results on stability of the bush B[2a] for the FPU-a model 
with those by Poggi and Ruffo obtained in for the FPU-/5 model [remember that 
B[2a] is the simplest bush consisting of only one mode i'N/2{t)]- The dynamical equation 
describing the evolution of the mode i'N/2{t) is the equation of the harmonic oscillator 
for the FPU-a chain and the Duffing equation for the FPU-/? chain, respectively. The 
problem of the loss of stability of the bush B\2a\ according to the parametric excitation of 
the other modes reduces in the linear approximation to studying the differential equations 



, [fl| ) for the FPU-a case and to studying the Lame equation for the FPU-/5 case. In 

both cases, coefficients of the above mentioned equations depend on j /N only, where j is 

a number of the excited mode and is the total number of particles in the FPU chain. 

Let us consider the main difference in stability properties of the bush B\2a\ for the two 

discussed mechanical systems. There are different threshold values of the energy localized 

in the mode Vn/2 for exciting the different modes Vj for the FPU-/? chain and, moreover, 

the modes for which sin^(^) < | cannot be excited parametrically at all 0. In contrast 

to this, we found that in the FPU-a model all modes with j 7^ ^ excited by the 
Note that the threshold of the loss of stability of the ZBM, i.e. of the bush B[2a], was obtained in 
Jl9[ and in the present paper with the aid of the Floquet theory. 
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mode z/jv/2 simultaneously, i.e., the same threshold value a Ac corresponds to them. When 
I aA |> aAc all modes ^ begin to increase at once, and, as a result, the original 
bush B[2a\ loses its stabihty. 

The studying of stability of the bushes other than B[2a\ for arbitrary N (which are 
multidimensional) is the more complicated problem, and we will not discuss it in the 
present paper. Let us point out only one interesting result. As was discussed in Sec. 
2, every bush possesses a certain symmetry group G and contains a number of modes 
whose symmetry is higher or equal to G. These modes Uj are depicted as sticks in 
Fig. 3a for the case of the bush B[8a] and we will refer to them, using the spectroscopic 
terminology, as "lines" (note that some modes of the bush are too small to be shown on 
this picture). The height of each such stick (line) is determined by the value of z/j(t) at the 
instant t. When N becomes sufficiently large (using the program PI, we considered the 
cases N — 128, 256, 512, 1024) a number of satellites, forbidden by the above mentioned 
symmetry restriction, appear near every line of a given bush (see Fig. 3b) for the amplitude 
of the root mode exceeding a certain critical value. This phenomenon may be described 
in the following manner. 

Let jo correspond to one of the modes uj^ belonging to a given bush. Then several 
"forbidden" modes Vj with j = jo ± 1, jo ± 2, jo ± 3, . . . arise above a certain threshold 
for aAc of a root mode. It must be stressed that the appearance of such "satellites" of 
different orders is not connected with any symmetry related principles. The excitation of 
the satellites brought about by the approach of frequencies of the appropriate modes to 
each other when N ^ oo and this phenomenon, in turn, results in resonance interactions 
between modes. The time evolution of the satellites can be very nontrivial, in particular 
for large N, amplitudes of satellites of higher orders can be greater than those of lesser 
orders. In Fig 3.b, we show schematically the satellite structure of the bush B[8a] for 
N=128 at t ~ 2500Ti, where Tl is the period of oscillations of the root mode. Note, that 
"first" modes, i.e. modes i/j with j near the beginning and near the end of the interval of 
permissible values of mode numbers, are also present on the Fig.3.b. The appearance of 
the satellite structure corresponds to the beginning of the process of the loss of stability 
of the given bush, but this structure can be observed for a sufficiently long time as a 
metastable state. 
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In connection with the above mentioned metastable states, let us note that there 
are a number of papers devoted to studying the dynamics of nonhnear chains above the 
threshold of the loss of stability of the zone-boundary mode, i.e. the bush -B[2a], in our 
terms (see, for example, and where some types of localized states were revealed 
in the certain cases). In particular, the very interesting phenomenon of the emergence of 



"chaotic breathers" as a kind of metastable dynamical regime was reported in ||2^. We 
do not study such phenomena in the present paper. 



4 Conclusion 

In the present paper we demonstrate how to obtain the invariant submanifolds in the 
modal space of an A^-particle chain with periodic boundary conditions by regular group- 
theoretical methods based on the concept of bushes of normal modes developed by us in 
[7-11] for arbitrary mechanical systems with discrete symmetry. Bushes of modes appear 
under certain initial conditions and they can be considered as geometrical objects with 
their own symmetry and as dynamical objects, as well. For the last case, differential 
equations corresponding to the dynamical systems whose dimensionality is frequently 
essentially less than can be obtained. The energy localized in the very definite set of 
modes, corresponding to a given bush, remains trapped. The object of such a type were 
revealed recently for the FPU-/5 chain by Poggi and Ruffo 0. We show that these objects 
do correspond to our concept of bushes of normal modes and that the group-theoretical 
methods developed by us earlier allow one to find easily all these objects for arbitrary 
chains (not for the FPU models only). We also compare the classification of the bushes of 
vibrational modes according with the different choice of the parent symmetry group Go 
(it can be the cyclic group of pure translations or the dihedral group). 

For a concrete chain — FPU-a model — we derive the integro-differential equation which 
describes the dynamics of this chain in the modal space in the limit — oo. In contrast 
to the KdV equation obtained by Zabusky and Kruskal for the FPU-a model in the 
continuum limit in the ordinary space, the above mentioned integro-differential equation 
is correct for all wave lengths, not only for long waves. 

We study the stability of all possible bushes of normal modes for the case N = 12. 
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For the bush i?[2a] to which the pattern of atomic displacements with translational sym- 
metry twice that of the equihbrium state corresponds, we found the boundary of the loss 
of stability in the limit —>■ oo. It was found that the bush B[2a\ loses its stability with 
respect to all other modes simultaneously when the amplitude of atomic displacements 
exceeds the value (0.303/a). Let us stress that this result corresponds only to the FPU-a 
model. Indeed, Poggi and Ruffo found essentially different results for the FPU-/3 chain 
(we give the appropriate comparison in Sec. 3.2 of the present paper). We also briefly 
discussed the appearance of the metastable satellite structure of symmetry forbidden lines 
near the permitted lines of the bush i?[ma] where m 7^ 2 for sufficiently large N. 
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Fig.l. Splitting of Eqs. (68) 
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Fig. 2. The monatomic chain in the equilibrium state and the pattern of atomic displacements for the bush J5[2a] 



(a) 



(b) 



Fig. 3. Satellites near the modes of the bush B[8a] (schematically), 
for = 32. (b) Appearance of the satellite structure for N = 128. 



(a) Absence of the satellite structure 
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Table 1. Irreducible representations of the cyclic group Go for N — 12 



Irreps 




a 


a2 


a3 
















all 


oymmetry 
groups Gk 


Modes 


To 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


[a] 




Ti 


1 


7 


7^ 


7^ 


1' 


1' 


7^ 


7^ 


7^ 


7^ 






[12a] 


^1 


T2 


1 


7^ 


7^ 


7^ 


1' 




1 


7^ 


1' 


7^ 


7^ 




[6a] 


^2 




1 


7^ 


7^ 


7^ 


1 






7^ 


1 


7^ 


7^ 


7^ 


[4a] 


^■i 


r4 


1 


7^ 


7^ 


1 


1' 


7^ 


1 


7^ 


7^ 


1 


7^ 


7^ 


[3a] 


1/4 




1 


7^ 


^10 


7^ 


7^ 


7 


7^ 




7^ 


7^ 


7^ 


7^ 


[12a] 




Te 


1 


7^ 


1 


1' 


1 


7^ 


1 


7^ 


1 


7^ 


1 


7^ 


[2a] 


V% 


Tt 


1 


1' 


7^ 


1' 


1' 






7^ 


7^ 


7^ 




1" 


[12a] 


V-j 


Tg 


1 


7^ 


7^ 


1 


7^ 


7^ 


1 


7^ 


7^ 


1 


7^ 


1' 


[3a] 


1^8 


Tg 


1 


7^ 


7^ 


r' 


1 


7^ 


7^ 


7^ 


1 


7^ 


7^ 


1' 


[4a] 




Tio 


1 




7^ 


1' 


1' 


7^ 


1 




7^ 


7^ 


7^ 


1' 


[6a] 




Til 


1 


7" 


^10 


1' 


1' 


7^ 


7^ 


7^ 


1' 


7^ 


7^ 


1' 


[12a] 





Here 7 = e^^'/^^ and, therefore, 7^ = i, 7^ = -1, 7^ = -i, 7^^ = 1 



Table 2. Bush classification in accordance with the parent groups Gq = T and Go — D 



Dim 




Gq — D (new bush) 


1 


[re] 


[14-18]; [rs-Tg] 


2 


[r4,r8] 


[r6,r3-r9]; [r2 - Tio, r4 - Tg] 


3 


[r6,r3,r9] 


[re, r2 — Fio, r4 — Fg] 


5 


[re, ^2, Tio, Fg] 


[Fi — Fii, F2 — Fio, F3 — Fg, F4 — Fg, F5 — F7] 


6 




[Fe, Fi — Fii, F2 — Fio, F3 — Fg, F4 — Fg, F5 — F7] 


11 


[ri,r2,...,rn] 
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